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\ 1. INTRODUCTION 


SOVEL was written as part of the Ocean Thermal Structure 
Analysis package for Fleet Weather Central, Rota, Spain. [Its 
primary purpose is to extract the temperature and salinity 
levels (by one-degree squares) in desired locations from an 
ocean climatology ee It takes the corresponding analyzed 
sea-surface temperature and the mixed layer depth at given 
locations which are read from input cards and adjusts the upper 
standard levels to these parameters. An additional feature of 
the program is the adjustment of the temperature and salinity 
gradients below the mixed layer depth. These gradients can 
be sharp indeed in the Mediterranean. Their sharpness varies 
with various parameters and seasons and this variation has been 
taken directly into consideration with this program. 

The program outputs are (a) the interpolated oceak (salin- 
ity and temperature) climatology (optional); (b) the tabulated 
adjusted ocean climatology (salinity and temperature) and 
sound speed (Figure 1); and (c) graphical printing of the above 
parameters (Figure 2). 

The subroutines MEDCLM and INTRP were written by Mr. Roger 


Bauer and the climatology tape was also prepared by him. The 


program will run on any CDC computer (CDC 1604, 3100, 6500, etc.) 


with FORTRAN IV or FORTRAN Extended compiler. 


40.0% LONG, 6.0€ 


TEMP, SALIN, VELOC, 


15,90 37.67 1513,50 
15,90 37.47 1513.98 
15,90 37,67 1514.47 
15,90 37.67 1514,680 
15.43 S7a7s 1513,59 
14,96 37,76 1512.35 
14,02 38,99 1510,13 
13,38 38.17 1508,55 
13,06 38,27 1508,44 
13,07 35,36 1509.41 
43,12 38,42 1510,47 
£3.17 33,46 1512.34 
13.14 38,45 1513,88 
13,07 38,45 1515.30 
13,01 38.43 1516.73 
12.96 38,43 1518.22 
12.95 38,42 1519.84 
12.94 38,42 1521.46 
12,94 38,40 1524.77 
12,96 3A, 40 1529.85 


Figure 1. Example of tabular printout of salinity, temperature 
and sound speed profile. 
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2. INPUT DATA AND REQUEST CARDS 


The climatological input data is available on a special 
tape prepared by Mr. Roger Bauer from data provided by Mrs. 
Margaret: Robinson of Scripps Institution of Oceanography. The 
tape includes data from the Mediterranean Sea, Black Sea, and 
the eastern Atlantic from the coast of Portugal, Gulf of Cadiz 
and Moroccan coast to 15° west longitude (20 to 50N; 15W to 45E). 

The synoptically analyzed sea-surface temperature and 
mixed layer depth data are read from cards in Subroutine J02. 
The first input card contains the number of profiles (LOP, 
Format I3). (See listing for Subroutine J02.) 

All sea-surface temperature data follow on one or several 
cards in Format 10F6.2 and the mixed layer depth is provided 
in the same format. (See listing for Subroutine J02.) 

The location cards (one card for each position) are read 
in Subroutine MEDCLM. The card format is: NLAT, NLATH, NLONG, 
NLONH, MONTH, IDAY. NLAT and NLONG are latitude and longitudes 
in whole degrees in Format I4. NLATH and NLONH are indicators 
of north/south latitudes or east/west longitudes (i.e., N, S, 

E or W only). The month and day are given with two numbers 


each in Format 212. (See listing for Subroutine’ MEDCLM.) 


i aa aan Wi a ls a deg eh 


3. 


ALAT 
ALONG 
DEDE 
DEP 


DEPTH 


DEPTH2 


EDP 


ICNT 


KE 


LINE 


ND 


LIST OF ESSENTIAL ABBREVIATIONS IN THE PROGRAM 


Latitude of the profile 

Longitude of the desired profile 

Depth in kilometers 

Depth of the levels in meters 

Depth of a given level in upper layers (above 
150 m) 

Depth of a given level in lower layers 

An additional depth counter in sorting of 

the values of a profile 

A counter for records per page 

A counter for the number of levels in a profile 
Number of characters in a line in printing a 
program 

A counter 

A counter for the number of profiles under 
computation 

The number of profiles requested by the program 
The number of values for lower layers for which 
annual mean values of salinity and temperature 
are available 

Indicator for north or south latitude 

Indicator east or west of longitude 


A counter 


VvoT 


vsti 


The number of levels in upper layers 

A counter 

Mixed layer depth in meters 

Salinity in 0/00 

Salinity at a given level in lower layers 
Salinity at a given level in upper layers 
An intermediate storage for depth values 

An intermediate storage for salinity values 
Sea surface temperature in degree Celcius 

An intermediate storage for temperature values 
Temperature in degrees Celcius 

Temperature at a given level in lower layers 
Temperature at a given level in upper layers 


Sound velocity meters per second 


An intermediate computation of sound speed 


4. SUMMARY OF ROUTINES BY FUNCTION 


] 
| 
| 
| 


SOVEL This is a control program which sets the counters 
LO (the number of profiles) and KE (the number of 
levels in any given profile). It calls several 
subroutines and checks at the end that all 
desired profiles have been computed. 

JO2 This subroutine reads from the cards the number 
of profiles required in the particular compu- 
tation, and the sea~surface temperature and the 
mixed layer depth for each profile. Calls 
Subroutine MEDCLM. 

JO3 This subroutine computes the sound speed for 
each level in each profile. The latest version 

of the sound speed computation formula from 

SACLANT (NATO) ASW Research Center, La Spezia, 

Italy is used. 

J04 This subroutine first prints, in tabular form, 
the month, day, latitude, longtude and the 

depth, temperature, salinity and sound speed 

for each profile. Thereafter, the values are 

graphically printed out for checking and 


eventual subjective correction. (See Figures 1 


ee ee 


and 2) 


a 


lieroy, Cc. C. Development of simple equations for accurate 


and more realistic calculation of the speed of sound in sea water. 
SACLANT ASW Res. Centr. Technical Report 128, 1968. 
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SOR This subroutine sorts the various levels in 


the profile after additional levels have been 


inserted below the mixed layer depth by Sub- 


routine GRAD 

GRAD This subroutine adjusts the temperature and 
salinity gradients below the mixed layer depth 
and adds additional levels, if so required. 

MEDCLM This subroutine extracts climatologicai data 
from a special ocean climatology tape. The 
tape is sorted from north to south and west to 
east. Data to be printed is selected on the 
basis of a request card (locations on the card 


must be in the same order as the tape). The 


a 


coordinate system used in sorting the tape is as 
follows: Longitudes are sorted starting with 

0° east and ending with 0° west. Latitudes are 
sorted with all data from a given latitude 


appearing before any data from the next southern 


latitude. 


5. LISTINGS OF PROGRAMS AND SUBROUTINES 


Listings for each of the programs and subroutines are 


presented in the same order in which they are summarized in 


section 4. 


PROGRAM SOVEL 


3200 FORTRAN (3,0)/RTS 05/25/73 


PROGRAM SOVEL 
NIMENSION TEMP(23),DEP(23),SAL(23),VEL(23),SS57(25),PL0(25),EDP(25, 
1O) LIVE (119) 2SDEP( 32), STEMP(32),SSAL (32) .NNN(25) 
COMMO. TEMP, DEPsSAL» VEL» SSTePLNSEDP.LINEs»SDEP»STEMP»SSAL,NNN,NUMSL 
TOR LE eb eRe SREP. NC MEAT HRT Hs SEONG: NLONHEASN THe IOAN 
LYal 
KEal 
3 call Joa 
CALL GRan 
CALL SOR CONTROL PROGRAM 
CALL JO3 i aaa a 
CALL 04 
LOsLOet 
IF (LOP-LO)5,3,3 
5 STOP 
END 


1 


SUBROUTINE J02 


© JO? 

(y 

( 

c 

¢ 

G 

(e 

C 

Cc 

C 
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: 

(> 
21 
23 
4p 
28 


3200 FORTRAN (3,0)/RTS 05/25/73 


SUPROUTINE Jd2 

DIMENSION TEYP(23),DEP(23),SAL(23)4VEL (73), SST( 25) ,PLD¢25) ,ELP (25, 
10).LINE(110).S0EP(32),STEMP(32),SS4L¢32, NNNC25) 

COMHON TEMP, DEP, SAL» VELsSSToPLIsEDP,LINEs SDEPsSTEMP» SSALsNNNg VUMeL 
LOP,LT,LO.VO,KE.KEP,NUGANLAT»NLATH,NLONG, NLONH, MONTH, IDAY 

READING OF VALUES 

MAX 20 PROFILES, 27 DEPTHS EACH, 

NUM IS NUMBER OF DEPTHS IN THE PROFILE 

LOP IS ‘IUHRER OF PROFILES 

POc IS PROFILE NAME 

nee nEPTH IN METERS 

TEMP IS TEMPERATURE IN DEGREES ¢ 

Sal IS SALINITY IN PROFILE 

ALAT IS LATITUDE OF THE PROFILE 

LTS IS TEMP, SCALE IN PLOYTING 

LTS 210 THEY O TO 10 DEG, IF LYSs4 THEN O TO 256 DEG, 

SST IS IN CENTIGRANE 

PLN IS THE MIXED LAYER DEPTH IN METERS 


FORMAT(I3) 

FORMAT(LOF4,2) 

1F¢L0~"4)49,40,28 

v031492,9 

READ 21, LIP : Reading the number of profiles 
READ 2%, (SST(]).121,L0P)> and SST and MLD for these 

READ 23,¢(PLD(1),121,L0P)) profiles. 

CONTINUE 

CALL MEDCL™ <«————- Calling subroutine for extraction 
RETURN of climatology from tape. 


END 


21 


26 


SUBROUTINE J03 


3200 FORTRAN (3,0)/RTS 05/25/73 


SURROUTINE 593 

MIMENSTON TEMP (23) .UEP(23).SAL(23),VEL(23),S5T(25).PLD(25),ELP (25, 
19),LINE(110),SDER(¢ 52),S5TEMP(32),SS4L (52) ,NNN(29) 

cOM4ON TEMP, DEP, SAL, VEL, SST,PLD, ED", LINE, SDEP,STEMP,SSAL, NNN, NUM,L 
LOP,LT,LILV9, KE, KEP. NUCANLATSNLATH, YLANG, NLONH, MONTH, IDAY 

c JO’ COMPUTATION 

ALAT = FLOATF(CNLAT#9)/10) COMPUTATION OF 

0 30 Il2*€,KeEP “SOUND SPEED 

IF (DEP (1)91,)500,900, 901 saa pepe 

Son veP(t) = 1, 
501 NENE =PEPC(l)/1000, 

VOT  =Se (TEPC L910, =O, ON Ge (aasreTeEMP({)-10, wets al, NGe(agse( 
LTEMP C1) 918, ) 072) 41, 2e(SAL(12 995, eM Adel CTEMP( 12918.) * (SAL C1) 255,) 
2y+NEP( 1) s61, 

VOP=9,1*ABSF (DEDE) **24(0.,0002*ABSF (DEDE)**2) *(ABSFCTEMP(1)=18,) 
Lee?ye( 0, JeneneraLat) 7/90, 

VOS 2, 0G = 7eTEMP( 1) #(ARSF (TEMP (1) 040,)**4) 

VSI =0.0N15#(ABSF(SAL(1)935.) 42) e(1°DEDE 

30 VEL(I)SVOeVOT «VIP +¢VOS «VSI 

ReTuRy 

END 


SUBROUTINE J04 


LOP,LT,LOeVIsK Ee KEP,NUCANLATsNLATHs NLONG,NLONH, MONTHS IDAY 
ALAT & FLOATFECNLAT®5)/10) 
ALONG © FLOATFCCNLONGRS)/40) 


9 3200 FORTRAN (3,0)/RTS 05/25/73 
SURROUTINE J94 70 
INTEGER XFIXF 
OIMENSION TEMP(23).QEP(25)+SAL C23), VELC23)eSST(27) PLDC 25), EUR I 25, 
16),LINF (410) ,SDEP( 32), STEMP (32) ,S5AL¢37),NNNC25) 
{ COMMON TEMP, DEP,SAL, VEL, SST»PLN ED? sLINE, SDEP.STEMPsSSAL NNN NUM SL 
} 


| f FO4 PRINTING OF ,VALUES 75 
42 FORMAT (25%, 5HDEPT41,6Xe54TEMP» OX, 4HSALIN, 6%, 6HVELOC./) 76 
43 FORMAT (25%,F5,0,6K,F 5.22 OX sF2,200%,F7,2) 79 
251 FORMATOD5X. 4HLAT . oF 7,40A2n4XsDHLONG, FP 7.1 ,A20/) 
253 FORMAT(1H1,25x%, 4HDATE, 16, 4%,19,777) 
E KK E SKE 
| PRINT 2535 10NTHs [DAY ) Printing of 
Bi PRINT 251, ALAT, NLATHs ALONGs NLO“H D, T, S, V in 
PRINT “2 : tabular 83 
PRINT 63, (08Pc 1), EMPcLySALel), VEL (1), LeKE,KEP)) form 84 
Cc GRAPHING OF VALUES 
59 FORMAT (tH4.29X%e75 GRAPHIC DISPLAY MF TEMPERATURE. SALINITY, AND 
TVELOCITY CHANGE WITH DEPTH e ///) 


261 FORMAT CL5x, tHT, 3K, 1H5, 198K, 2410,24X,2H15,18X, 25120, 18X,2H25,18x, 2430 


) 
a2" ¢qnMar ttAx sine 2e, Susie tay, cae, 204) 8nds Dens Ons, tax, 2uSa, lex, 
1240) j 
53 FORMATC19X,1HV,35K,1540, 18%, 2H20,18X%,2H499,18X%,2HO0,19xX,2HB0,17x, 
LSH409) 
PRINT 56 ~, 
PRINT 264 
{ PRINT 52 
S| PRINT 53 
{ ITSLHT 
a | PSelus 
IVetHy 
1X2LHx 
lyeLHy 
IBLANK=41 
[Polis 
[Melle 
T7217 
DO 64 Is1°100 
64 LINEC])=SIM 
no 66 1210,100,10 
66 LINECT)e#]P 
PRINT LAa%,CLINECT)2fe1.100)9 


11% FORMAT(19X,444,100A1) Graphing of 

DO 68 121,99 Diy Ty SV 
6R LINECT)#IBLANK by printer 
74 1eKE 


TF (DFP (1 )=150,)75,75,72 
72 VEVUISDEPCI) DEP (1-4) 
IF CDEP CT) 91500,)750,750, 720 
j 720 KusxXFIXFCUFUTs100,) 
6u TO 74f 
| 73n KLaXFIXFCDEUL/50,) 
) 740 KaAKz1 
73 PRINT 200, ¢(LINE(L)»Le1.99) 
200 FORMATS LR pa eg ggeesne) 
IF CKAK KL) 74,72, 79 
74 WAKE AK ed 


SUBROUTINE J04 (continued) 


ko TO 73 
TEMPERATURE 
90 210 t=10,190,19 
LINE(L)=IP 
JeXF PXF (TEMP( 1) 04,2220 
1F (5982.82.83 7 
Jat 
IF (J~w100)85,84,84 
J299 
LINECJ) SIT 
SALINITY 
eXF XE (SaL C920, .2 500) 
1F0M)92,92,9 
Met 
IF €%e100)95,94,94 
Ms,99 
TF (M#J197 296097 
LINE CAI SIX 
Go TO 100 
LINE(M) 215 
VELOCITY oe 
TORRE TEE GG (1) 9145 Graphing 
1F (Ju 102,102,105 by printer 
JJel 
IF (JJe100)105,104,104 
JJR99 
Te (aden 20a0206) 198 
LINECSJETY 
Go To 119 
LINE(JJ)3 Iv 
IF (Joys) 124,112,114 
Tae 
IF (DEP C1) ) 400,400,124 
DEPCI)22,5 
IDEPSXEIXF CDEP CT?) 
PRINT 310, 1DEPe(LIVECT), 181.99) 
FORMAT (10X%,18,1K,14e,99AL, LHe) 
NO 333 1=1.99 
LINE ])s1BLANK 
KEekKEel 
IF (KE @KEP) 71,71, 345 
KESKKE 
RETURN 
END 


1 


ee 


SUBROUTINE SOR 


3200 FORTRAN (3,0)/RTS 


SURRMUTIMFE SOR 


C572°/73 


NYMENSTOM TEMP C23). DEP (23) .SAL 625). VEL 623) eS5T( 25) ePLU! 25) EL PC 5, 


1°), LINE (119),SDEPC $2), STEMP C32), SS4L (32) NNN 22) 
ST,PLO,EDP, LINE, SUEP, STEMP, S5ay ,NNN,NUM,L 


“OMMON: TEMP, DEP, SAL, VEL, 


LOP.LLT,LI,VI.KEDKEPSNUCINLATANLATHs YLONG, NLOMH,MONTH, I DAY 


LuUsd 
TF CEDP CL, 29100, 1504.4 
INN(LUIANU 1e2 


90 '0 8 SORTING OF T, S VALUES IN 


IF CEDP(LO,5))50506 

INN(LO) aNU4eS 

Go To 8 

re CNUC#1917,16,17 

INV (LO) SNUM+3 

30 TN 8 

MNN CLO) SNUM eG 

KKFSKE 

NMNSKKES+SNN( OQ) ad 

KEPaNin 

LLMSNHN(LO) 

po 295 yot,tum 

IF ¢L-1)11.11.12 

SpePcpy=n,1 

STEMP(1)=SST(LO) 

KEzKee! 

RQ TO 105 

[FeLi 13914,49,13 

TF éLUH-3)40.,70,40 

IF (DEP CKE) @=PLD(LO)) 19,241,20 
IF (CDEP CKE #4) »PLD(LO) 924,235.22 
SSALCI*4 2 S5AL (KE~1) 
SpeEPC1LI=PLA(LU) 
STEMPCT)=SSTELO) 
SSalL¢1)sSal(<eety 

Luss 

Go To 65 

DAL CLL) SCSAL (KEL) @SALCKEDI/2, 
SSAL CL) =SSAL CTR) 
SpEPCLY SPL (LO) 
STFMPC(IT)®SST(LO) 

Lusk 

KEnkt +1 

GO TM 19% 
SSALCLeq) ¢SAL (KEWL) eSALCKEII/2, 
SSALCL=ESSALCTeL) 
STEMP (1) 2S5T (10) 
SVEPCI)I=SDEPCT) 
SpePCr+t)=PL (Ld) 
STEMP C141) =STEMPEL) 
SSALCLo1)25354L¢1) 

KEshbel 

Jeet 

Lus 

gu to 105 

SSaLCpelrye Sai (KE=l) +SAL (KE) SAL CKE S15, 
SGALCLYSSSALCT #4) 

SSALCL+1)s59aLe1) 

STEMP CT) SS9TOLO) 

StemP¢ tel) eSstcLa) 


SHEP CL) Spee (Ke) 


CORRECT DEPTH SEQUENCE 


pagiti min Apes wt clag sim 


SUBROUTINE SOR (continued?) 


SOFP (141) 2PL9(L9) 
Tel+1 
KEsKEe? 
Lyset 
Go To 105 
24 SSALC LoL) 2CSAL (KE el) @SAL (KF) SAL (KF 41) 9/3, 
SSALCL) sSSAL (fed) 
SSalclertd#SSaccl) 
STFEMPC(TIESSTI(LO) 
STEMPC1e1)2S5T(L9) 
STEMP C14?) 2S55T(LO) 
SSAL (1 *2)25SaL 01) 
SpFPt]) epEP (Ke) 
SpEP Cel) anEP (Keel) 
SpEPCL+2)aPL (LO) . 
1F (DEP {KE*2) -PLD(LO))28, 76,29 
2° KESKEe2 
RO To 27 
26 KESKEs3 
27 12142 
Lust 
30 To 105 
28 SDFP (142) = EP (KES?) 
STEMP( 145) 2S58T(LO) 
SSaL(1*3)=SSaL C1) 
SDEP(] 9) sPLU(LO) 
IF (CDEP CKE*3) «pLD(L9)) 52,50, 29 


29 KESKFe3 
GO TO 31 

3) KESKE+4 

Si fste3d 
Lust 
GY To 105 


32 SSaLCreLy=(SaL (Keel) SAL (KE @SAL KE 4!) Sal (KES2 eSab (KE SS)I/5, 
SSALCL) SSSALC Let) 
SSalL¢l+ly=ssal (1) 
SSALCL #2) 2S835aL (1) 

SSALCL*#5) 2SSaL (1) 
SSALCL*#4)S5SaL (1) 

STEMP (144) 2SST(LO) 

SDEP (1 +3)snEP(KEs9) 

IF CDEP (KE +4) @PLO(L0)) 56, 54, 53 

53 spepc1+4)s°Ln(L9O) 

KEsKEs 
ud TO 35 

34 SPFP CI +4) SHEP (KE+4) 

KEskEe5 

35 rates 

[ LUs 

$6 SPEPCI+4) =nEPCKES4) 
SSaL(1*5)esSaLite4) 
STEMP( 145) 23ST(LQ) 
IF EPLOILO}< 2004} 62,3561, 361 

361 PLO(LI, S209, 

362 SUFPCT+5)2PLD(LO) 

LF COEP (KE SS) ePLD(L1))38, 58,57 
37 vEskee 

30 TN 39 
3A KESKE S46 
39 [ste5 

Lust 


— sete teen 


{ SUBROUTINE SOR 


4? 
45 


| 
4} 
; 


(continued) 


105 

TE CENPILO,2))700700 41 

IF FENPCLO,2)—100,)42, 4%, 44 
1s CENP(LO,5))463,43,45 
SUFPCTISEDP(LO.1L) 
SDEPC[ +1) eepP(La, 3) 
STewPrtyeknP¢, 9,2) 

STRFMP C11) sENP(L 9,4) 
SDISCSSAL (Te1)*SALCKE) I<, 
SSaccpye5Sal(tel)«Sp! 
SSALCI*#E) SSSAL Clad ie 2, eSif 
Toto 

LUsg 

OG TO 405 

SDFP C1) SFDP(LOs1) 
SUrPCi#t)SEDP(LO, 5) 
SDpEPCl+2)sEDP(LO,5) 
STEMPCIYFENP(LO.2) e400, 
STPME (ied SBOP EL Ds”) 

STEM? (142) sEDP(L0,6) 
SDIS(SSALC1=4)2=SALCKED) 24, 

SSALCTI SSS abd | sleSpr 

(let) sSSAal( lei .e2, 051 
1*2)S5SAL CJ o o1)83,eS)] 


10 T 


SSat 
SSAL ( 
fe. 

WEsKF +4 
tues 
Go Ti 105 
SDFP(I)SEDP(LOs4) 
SUEPCI+1)2EDP(L0,9) 

Scene. +2)2EDP (10,5) 

MPCTSFENP(LO9D) 
Srrieet iene 6 es 
STEMP (142) sE0P(L0,6) 
Spa CSSAL CIT 1)@SALCKEDI/4, 
LOT ESSaAL eC] @t)9S01 
SSaL clei) eSSar (let) o2, eS] 

TSAL CI #2). 55SAL (fet) e4, 259] 

[a]? 

LUs3 

66 TN 105 
SDFPC]) apeo KE) 

Sahel dod, ner TemMPc(wpy) Stn, ib a 
STEMP el) STEMP(Laided, as 
ho TN $42 
STEMPCT)STE4P (KE) 

TF tSSAL (101) ~0,25~9au (er) 324, 344, *1% 
VF t$520344,314,402 
songs! LLM} 492¢4g2¢334 

LOPPSCSGAL CL #2. eGaL tke?) 972, 

4. i 315 
eau 't sSar 
KEKE #] 
COMTT JUE 
KESKKE 
Ll 
ju 291 peKe, My 
TEMP CT SSTEMP CL) 
DEP CT eS0EP CL) 
SAL CTI *5SAL CL) 
Lest. 

CUT yp yur 

RETURN 

END 
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(ke? 


SUBROUTINE GRAD 


a a 


NPCs >= 


30 TM 14° 
tc ne Tweet LF 
TE (PL UC) 


Qms » 
1OP,LT, m , 
!zl) 
UCse 
( 1% fare 
TFC UEP CL r=, 
SA NUCS1 
B IF ENeP eyed) 
41 ius 
6 TF CORP CT ye 00 
~* TR Cpe +} 
nO ce ‘ 
t Fe* a | 
ames ( jy#9 
7 wr pee 
, e035 Way 
t SST ) 
’ me 7 ss 
> nm? F 
o “ . 
3% ifs 
Oo ‘ BOG « 
31 PSC FRE 
(J,2)2587T 
j ty T 
VP CI,9): 
3U 115 
9 IF SST J) 
4 . 
IF CDIF*O.499 
gn f SAl 


» +0 
a3 
ty) 
gi 
ta) 
(gy) 
(J) 


T2N0 FAxTRAN CEL LRTS SA 2° £73 
PEN DEPE ZS) pGAL C26 CEL FT 2B. ESET QR ISPLO! 292 6b VOC 6, 
£°( 39), T MPC S2Y,S02L S22 .NNYE27) 

SAL, VEL, 28°, Pk", Fue, ve, SEs, STeMP,SSay NIN, UL 
KEP, JC Neat Nia Tw ING, “it £ i, “*INTHeIGay | 
ya }, 125 ADJUSTMENT UF GRADIENTS 

BELOW THERMOCLINE 
CI2F17V5. 42, 4 
(J) )G00"9 
{ big ae de te 

LOw aN? 40% SMace rSTawre 
’ »? 

Pepe 2) I74, 

13383, 42.82 
{ 1)) 

HUF 

2. °DDF 

5,°DD 
efiFeggn 
“2,¢DIF 
e5,eDIF 
PC Tags 
0,94 
(1 93.92.92 
Lat sypv*, 
eONF 
#2, e)0F 
+3, enor 
1?1 

i) 
w2,e119, 

e) DIF 
alts too 
©2.°DIF 

y 5 Dt 
a) AND 1M SMaue LeTance 
3226426 


pOPY AVAILAPLE TO DDG DOES """7 
PEnciT FULLY LEGIBLE PRODUC: J. 


SUBROUTINE GRAD (continued) 


DLFSCSSTCS)=TEMP CT o1)974. 
} TF CDIFSO,2927012700271 
27 TF CSAL CT =1)"SAL01)2275,272,272 
272 N[F=n, 
‘ GO TO 271 
273 ULF zg. 0% 
271 EDP J,2)=8STC Sent Fel, eSeylEr 
EUP CJ. 4)SENP( Je 2 en Fel, lenir 
FDP(J,6)=2E0P(J,4)00,8eDI© 
G0. TO 415 
| 20 NDFSC IEP (Ley ePLNC 19746, 
FDP J,1)2PLN( Uy) +npF 
EDP CJ,3)2PLD (J) 42, eDDF 
EDP (J*5) 2PLD(U) +3, DNF 
DIFSCSSTCU)eTEMPCTe1))74, 
IF (D1 +*0,1)290,280,281 
280 LF CS4 C194)" SAL(61)928392%909282 
2g¢ olr=0, 
go 1 284 
3 NUFsey,0A 
1 EpPCU,edESST Cy) =nI Fe, Senter 
EDP (de 4) 2ENP (de 2d aN Fed, lenir 
EDP (J,6)=E IP (J, 4) oD, PODIF 


60 TO 115 
| c MLD APPROX, 3ETWEEN TW LEVELS, 
35 NDFSCJEP (Led) -PLNCIIIZ3. 
H EDPCU, 1) ®PLOCJ) epoF 
FDP (5,3) =PL0(N) 2, DOF 


EDP (J95)20, 
| NEFECSST OJ) TEMP( Jal) ) 7/3 
EDP(Js2)=SSTC sent Fed, eure 
Fp? (4,4) ENP ey, 2) -01F 

1415 CONTINUE 
RETURY 
END 


| COPY A To DOES NOT 
Hr iuthe ieuteul ~ S3R0SCTION 


ij SUBROUTINE MEDCLM 


| 3200 FORTRAN (%,0)/9TS 05/25/73 
sIROUTINE MEDCL™ 
; ¢ PROGRA™ MENCLM RErleves CLIMATOLOGY NaATA FRIY TAPE, 

G TAPE [S SORTED ORTH TO SOUTHe WEST YO EAST» 


ATA TO FE PAINTED I5 SELECTED ON FaSIS IF REQUEST CARD, 
HEQUEFST CAQDS MUST BE IN SAME ORNE? 45 TAPE, 


C THE COORD VATE SYSTEM JSED IN SORTING THF TaPE [S 4S FOLLOWS 
) ; LONGITUDES AWE SURTEN WEST TO EAST STARTING WITH OF AND ENDING WITH OW 
ia c LATITUDES aR SORTED NORTH To SOUT~ WITH ALL NaTa FOR pa LATITIInE 
FY Cc APPEAPING YFFORE ayy pata FOR THR NEXT MORE SJUTHERN LATITUDE 
| NEMENSE IN TEMP (23), DEPC23).5AL(29),VEL(29)4SS7T(25).PL0(25) EVP C25, 
| 19) LINE ( 119), SNEP(32),STEMP(32),554L032) .NNN( 25) 
|) NOMMON TE M2 NEP, SAL, VEL ,5ST,P L",6D°,) INE oF P,' Ty MP,S5Sac,NN Ye NUMAL 
' LOPLLT LO VI,eKE KEP, NUCANLAT,NL ATH, YLANG, VLONH, MONTHS IDAY 
5 
‘ 


QEMENSION NEPTH(O), TEMPS(12.5),S5aLTS¢ 4) 
PIMENSTAt YEPH2¢17), TEMPNC17), SAL T9497) 


C INTTHALIL72 COUNTERS 
| C TONT 15 cOINTER FOR RECORDS PER Pate, 
; NLAT,NLAT ARE Toe REQUESTED LaTITUpe sESIGNATORS 
c SET LOGICAL INITS FOR 1/79, IRFCARD REANER, [Ts TAPE BNET. 


JLATs JL .'Ged 5 ]T2i 
HEAD A QF IEST CAR 
READ 4907, UL ATs ILATHs NLONG sa NLONY, MONTH, [DAY 
1900 FORMATCI4,A4141,14,41,212) 
C SHECK FOR TYE LAST CARD 
IFONLATH.E ).iHn 2 STOP 
CHECK LaTITne Au‘upD LONGITUDE 
q MLATeNLa’ 
| IFENUATY, © ),LHND22, 
4 21 MLATSoMLAT 
e? is ell del AT 
MLONGENLONS 8 TP ONLONK, E92, 14W)23, 28 
23 MLONGHS oS ORAL UNG q 
25 LF e'4l g KAT AO, aq, tr 
3n be RESP a dL ats cate st 


27 


rc CARD OUR OF hadi > PRIVT ERROR MESSAGE, READ ‘Ht XT CARD, 
4n eae LOCL, nthe. N pATHe LONG MUON HNN Ty 1D 4 
1001 FORMaT( 2294 wn 9 F SEQUENCE, 3X,14,41,14,41,212) 
} AETURN 
C READ NEW TaPE RECORD, AND CHECK aGatnNST REQUEST card 
ee 
34 REAUCT eLO92)LATALATHs LONGe LON, S SaND 
1¢nF tetonctenPSeys ta, Ja1,12),7ALT% (1), 1 et, 9S) 
STUER ECL? FEMPDE LY gSaLTOCl irr 
4 1902 FOoRMATC?A¢T4, 2410021975 (F8, 2.13F4 ere eet ae 
wackSpach 1° 
NUM @ NS @ YD é 
KEP = NIM 
e 
CHECK TaPE RECORD LaT aN LONG, 
JLAT = LAT B LFCLATH, EQ, 1HNI52,52 
51 JLATSeJLAT 
52 sratsgngeuLal 
JILONG = LONG 
IF CLO, bo, law 55, 2 
55 JLONT 2 3600~ jLONG $8 a0 To 25 
rc NEFD Tape STATUS CHECK ON ABNVE TaPETso, 


COPY AVAIAMIE TQ ppe DOES NOT 


Péci.iT FULLY LEGICLE P 


SUBROUTINE MEDCLM (continued) 


109 


400 


410 


411 


500 


210 


‘avE PROPER TAPE RECOR)!, 


adic erly REQUEST 


Cath INTHPCIONTH,LOAY, TE1PS. TEMPN,' SND, TEMP) 


190 400 J#1,NS 

DEP CT) SDEPTHCL) 
SALCL)SSALTS(1) 

IF ONS, GE,6)999 TO 414 

"NS es wS #1 

190 410 [=2NS,23 

pert reseic: sTEMP(])ae0, 
ETURY 

NDCSN0%6 

no 590 p27, NDC 

DEP CTY SDEPH2(1 9) 

Sac (1)=SALTD¢1-6) 
IDPCSsnpCel 

po 540 18 Dc, 25 

DEPC I) SSALCI) STEMPC(] 290, 
ENN 


RGUCTION 


al Ld gai Sr 8 VEEN AL rth) al lglg 


SUBROUTINE INTRP 


3200 FORTRAN (4,0) /RTS 25/28/73 


SURRHUTINE INTRP(MINTH,IDAY, TEMPS, TEMPN,NS,ND,T) 
NIMENSION TE 49017), TEMPS(12.9)4T0239,4(5),805) PAMTM( 12) 
1,09012,5),51012,5) 
DATACDAMTHS3S1 6 2Be a Slee 50,0 Lee Gee Ste Set S0 es SL ca Bead1.)s 
1¢RAD5O=.%233987755 4, (cot1)30, ) 


} c COMPUTF CONSTANTS JSED IN HARMONIC aNAL YSIS 
IFCCOCLINELO.) 150590 ‘ 
90 NO 190 [=1,12 $ 90 100 Kel,s? 


ANGSFLOATE (Ley) #RANSO@FLIATF IK) 
cCOCE Ky scOSF (ad) 
100 SI€1,K) S57 °F CANG) 


c FIND ANGLE Iv a dUAL CYCLE FOR THE GIVEN Day 
150 ANGSCCFELOATE ( [DAV @15) @DANTH(MONTH) Y/N AMTH (MONTH) @FLOATF (MONTH) Do 
1Ranso 
c cOMPUTE THE 12 4ARMONIC TERMS AND THEN FIND THe TEMP ON THE GIVEN Tay 
vO 400 121,'S 
AUsA%20, 


pO 469 K=1,5 
ACK)=B(KI=N 
no 149 taka? 
A(K)=A(K)#TEMPS (1, e CIC, K) 
160 8(K)=2B0K)*TEMPS( 1.) #SIECIT,K) 
po 200 t=1,12 
AQZANSTEMPS(T, 1) 
16 (172*2,E9,1)170,180 
179 AbzA%eTEMPS(], 2) 
60 To 20n ; 
180 AG=A64TEMPS(T,J) 
| 209 CONTINUE a 
TS) e(AdeAnwCO5F(AIG*6)) 712, 
UO 240 121.5 
ANZSANGHFLOATE (1) 
240 TESDETC UI *CACT) MCOSF CAN) BCT) #SINF CAND IZS, 
IF ¢J,E9,1)400,241 
241 IFCTOCJ),GT.T(521)) 245,409 
245 senontH 
IF (I DAY,6 7 (15)260,250 
250 M=Meo1 
1F(M, GT. 9) 769, 755 
795 4212 
260 IF e TEMPS (NM, 4), 37, TEMPS (M, dod) 4005779 
270 MeMey 
| _ IF OM,LT. 129279,275 
a?> vst 
278 IF ¢TEAPS(M, J) GT, TEMPSEM, Jw) ) 4004750 
280 TCJ)2TCUe1) 
400 conty juE 
po 4¢y pel, in ' 
429 Teles) 2TEMPOCT) 
vO 500 121,10 4 
te T6149) G7, 70%) 450, 40 P 
450 T(l+4)=T(6) 
e CHANGED 1 T9 [#6 [1 420+ 470420450 
509 CONTI WUE 
519  EeTURY 


| ENT 


